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Abstract. We present a systematic study of the effect of metallicity on the stellar spectral energy distribution (SED) of O main 
sequence (dwarf) stars, focussing on the hydrogen and helium ionizing continua, and on the optical and near-IR lines used 
for spectral classification. The spectra are based on non-LTE line blanketed atmosphere models with stellar winds calculated 
using the cmfgen code of Hillier & Miller 1 1998). We draw the following conclusions. First, we find that the total number of 
Lyman photons emitted is almost independent of line blanketing effects and metallicity for a given effective temperature. This 
is because the flux that is blocked by the forest of metal lines at /t < 600 A is redistributed mainly within the Lyman continuum. 
Second, the spectral type, as defined by the ratio of the equivalent widths of Hei A4471 and Hen /14542, is shown to depend 
noticeably on the microturbulent velocity in the atmosphere, on metallicity and, within the luminosity class of dwarfs, on 
gra vity. Third, we confirm the decrease in Teff for a given spectral type due to the inclusion of line blanketing recently found by 
e.g. lMartins etalj f20Q2). Finally, we find that the SED below ~ 450 A is highly dependent on metallicity. This is reflected in the 
behaviour of nebular fine-structure line ratios such as [Ne ni]/[Ne ii] 15.5/12.8 and [Ar iii]/[Ar ii] 9.0/7.0 yum. This dependence 
complicates the use of these nebular ratios as diagnostic tools for the effective temperature determination of the ionizing stars 
in H n regions and for age dating of starburst regions in galaxies. 

Key words. Stars: atmospheres - Stars: early-type - Stars: fundamental parameters - Stars: abundances - HIT regions - planetary 
nebulae: general 



1. Introduction 

Massive stars are a dominant force in the evolution of the in- 
terstellar medium of galaxies. The extreme ultraviolet photons 
of massive stars can ionize surrounding atomic hydrogen gas. 
This ionized gas is heated by the photo electrons and cooled 
through forbidden line radiation of trace species such as oxy- 
gen, resulting in a temperature of some lO"* K. Besides this en- 
ergetic coupling through the radiation field, massive stars can 
also influence their surroundings dynamically. The high pres- 
sure of the ionized gas will drive shock waves in their suiTound- 
ings, sweeping up the ambient molecular gas. Moreover, these 
stars also have strong stellar winds and explode as type ii su- 
pernovae, both of which provide kinetic energy to the interstel- 
lar gas. Understanding the characteristics of massive stars and 
their interaction with their environment is therefore a key prob- 
lem in astrophysics. 

The stellar spectral energy distribution in the extreme ul- 
traviolet (EUV) controls the ionization structure of Hii re- 
gions. However, because of the high opacity of neutral gas in 



the EUV, this wavelength range cannot be observed directly. 



Spectral typing of stars is generally done through optical (( 


5-g- 

)or 


IConti & Alschuleill97lVConti & Frost'l977VMathvs'l988l 


near-IR features (e.g. 


.Hanson et al. 


1996; Meyer et al. 


,19981 


iHanson etalJll998l 


Lenorzer et alJ 


|2002|). Alternatively, ob- 
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served nebular ionization structures can be used to probe the 
ionizing fluxes of the stars energizing the medium. This tech- 
nique is now widely used to determine the properties of stars in 
the EUV and has even evolved to a more general astronomical 
tool, paiticularly for the study of regions where individual stars 
cannot be directly tyt^ed (e.g.lOev et al ]l200(ilTakahashi etalJ 
'2000'; 'Ok amoto et ai1l200lt iMorisset e t al."2002'). Moreover, 
because the hardness of the EUV radiation field is a good 
measure of the spectral type of the ionizing star and because 
later spectral types live longer, the observed ionization struc- 
ture can be used as an age indicator of a starburst region (e. 
ICrowther et alJl999tlTTiornley et alJ2000llSpoon et alJ200i 

In recent years, some problems with the latter technique 
have emerged. For example, for the well-studied Hit region 
G29.96-0.02, the spectral type derived fr om direct observa- 
tions of the stellar lines in the ne ar-IR dWatson & HansonI 
ll997HMar tfn' -Hernandez et al. '2003^) indicates a much earlier 
spectral type than that obtained from the ionization structure 
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jMorisset et alJ2002i) . This seems to be a general problem and 
may be related to the metallicity of the region. Indeed, there is a 
loose correlation between the ionization structure - as measured 
by for example the [Neiii]/[Neii] 15.5/12.8 iim line r atio - 
and the metallicity of H ii regions (see Martin-Hernande z et alJ 
l2002hl) . Such a correlation may reflect the effects of metallicity 
on line blanketing or on the characteristics of the stellar wind. 

Much theoretical effort has been dedicated to best describe 
the EUV spectra of massive stars. This is a formidable task be- 
cause these stars have strong winds and extended atmospheres. 
This leads to strong non-LTE effects in the formation of spec- 
tral lines. Over the last ten years much progress has been made 
and current models include the effects of tens of thousands of 
lines on the energy balance and temperature structure of the 
stellar photosphere and wind. So far relatively modest eff'ort 
has been investigated in systematic studies of the effects of 
metallicity on the stellar spectral energy distribution. In par- 
ticular, there is no good theoretical understanding of the effects 
of metallicity on the ionizing fluxes of massive stars or on the 
optical and near-IR spectral characteristics used to type these 
stars. Here, we study the influence of metallicity on the spec- 
tral energy distribution of O stars and determine its influence 
on the resulting ionization structure of H ii regions. 

This paper is organized as follows. Sect. |2l presents a set 
of main-sequen ce (dwarf) star models constructed using the 
CMFGEN code bv lHiUier & Millet! ( 1 1998 1) and compares the pre- 
dicted EUV fluxes with those from other codes. Sect.|3]presents 
a detailed analysis of the variations of the EUV spectral appare- 
ance and ionizing fluxes with effective temperature and metal- 
licity. Sect.l^investigates the influence of metallicity and other 
stellar parameters on the optical and near-IR spectral calibra- 
tion. In Sect. |5j the ionizing structure of single star H ii re- 
gions is studied. Finally, Sect. |6l presents the conclusions of 
this study. 



2. Description of the models 

The models present ed in this paper a re con structed using the 
CMFGEN program ofl Hinier&Milleil ( Il998l) . to which we re- 
fer for a full description of techniques. In short: cmfgen itera- 
tively solves the equations of radiative transfer subject to the 
constraints of statistical and radiative equilibrium, for an at- 
mosphere with an outflowing stellar wind. The ions included 
in the non-LTE calculation are H i. He i-ii, C iv, N in-v, O iv-vi. 
Si IV and Fe iv-vii. The cumulative effect of the iron lines, caus- 
ing a line blanketing effect, is self-consistently accounted for 
by employing a description of the atomic models in terms of 
super levels. All and all, approximately 9000 transitions are 
accounted for 

We used this code to calculate a grid consisting of eight 
main-sequence (dwarf) star models ranging in effective temper- 
ature from 33 000 K up to 51 000 K. Basic stellar parameters, 
Teff, Rt and M,, using masses predicted by evolution theory, 
are taken from the Vacca et al. (1996) calibration and are given 
in TablelU 



We have opted to describe the photospheric density struc- 
ture using a constant scaleheight 



H 



(1) 



where yu is the mean molecular weight in atomic mass units, geff 
is the effective gravity at the stellar surface corrected for radi- 
ation pressure by electron scattering and all other parameters 
have their usual meaning. Near and beyond the sonic point, the 
density is set by the velocity law through the equation of mass- 
continuity. The connection between photosphere and wind is 
smooth. The wind velocity structure is given by a standard 
law adopting P = 0.8, which is r epresentative for dwarf O stars 
dOroeneweeen & LamersI 19891) . The terminal flo w velocity Voo 
follows from a scaling with the escape velocity ("La mers et alJ 
Il995.) and is al so listed in TableG] Themass-loss rates are from 
predictions by IVink etalJ (l200n[ 12001 h . These authors show 
that for O-type stars the mass loss scales with the metal con- 
tent Z as 

M cc Z°-^^ . (2) 

iLeitherer et al.l ll 19921) showed that the final wind velocity also 
depends on metallicity as v^o oc z" . We do not incorporate this 
dependence in the models. The microturbulent velocity was set 
to 20 km s"' for the atmospheric structure calculations. For the 
calculation of the emergent spectrum we assumed microturbu- 
lent velocities of 20, 10 and 5 km s"' . 

In order to investigate the spectral appearance as a function 
of metallicity, we have calculated an additional set of models 
for our models #2 and #5, in which we set the metal content to 
Z = 2, 1/2, 1/5 and 1/10 Zq. In conjunction with Z, we adjusted 
the heUum abundance Y according to 



^ Z. 
AZ 



(3) 



where Fp = 0.24 is the primordial helium abundance (* Audouzel 
[1987) and (AF/AZ) = 3 is an observed constant (Pagel 19_9_2) . 
To arrive at the a ppropriate m etal abundances, we simply scaled 
the solar values (Cox 2000), which we give in Table|2for ref- 
erence. The mass-loss rates and scaleheights, which also de- 
pend on Z through the mean molecular weight, of the additional 
models are given in Tabled 

2.1. The assumption of a constant scaleheight 

The assumption of a constant scaleheight could lead to an over- 
estimation of the efl'ective gravity as derived from density sen- 
sitive lines. To test the validity of this assumption we recalcu- 
lated model #6 usin g an atmospheric structure calculated with 
the ISA-WD4D code of lde Koter et alJ([l993lll99'7l) . This provides 
a more realistic density stratification as it accounts for the ef- 
fects of continuum radiation pressure on the atmosphere (see 
section 12. 2> . For the comparison we focused on the spectral 
signatures investigated in this paper, i.e. the EUV, the optical 
He I /14471 and He ii /14542 lines and the photospheric lines in 
the near-IR K-band. 

Close inspection of the predicted EUV spectra shows that 
the number of ionizing photons below a certain wavelength, as 
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Model 




M, 

'«Q 


Rt 

"o 


(cm s~") 




(10'*R^) 


logM 




(kms"') 


log Qo log Qi 
(photons s"') 


1 


51230 


87.6 


13.2 


4.139 


6.032 


7.703 


-5.209 


0.8 


3450 


49.87 


49.38 


2 


48670 


68.9 


12.3 


4.096 


5.882 


8.308 


-5.375 


0.8 


3240 


49.69 


49.18 


3 


46120 


56.6 


11.4 


4.077 


5.772 


8.375 


-5.599 


0.8 


3140 


49.50 


48.95 


4 


43560 


45.2 


10.7 


4.034 


5.568 


8.971 


-5.805 


0.8 


2950 


49.31 


48.72 


5 


41010 


37.7 


10.0 


4.014 


5.404 


9.046 


-6.072 


0.8 


2850 


49.08 


48.44 


6 


38450 


30.8 


9.3 


3.989 


5.229 


9.314 


-6.369 


0.8 


2720 


48.82 


48.05 


7 


35900 


25.4 


8.8 


3.954 


5.062 


9.712 


-6.674 


0.8 


2570 


48.51 


47.43 


8 


33340 


21.2 


8.3 


3.926 


4.883 


9.934 


-7.038 


0.8 


2450 


48.06 


46.07 



Table 2. Chemical composition adopted for the models with 
Z^Zq. 



Element 


Mass fraction 


H 


0.7023 


He 


0.2820 


C 


3.050 X 10-3 


N 


1.100 X 10-3 





9.540 X 10-3 


Si 


6.990 X 10-'' 


Fe 


1.360 X 10-3 



Table 3. Variation of the scaleheight and the mass loss rate with 
metallicity for models # 2 and #5. 





Model #2 


Model #5 






logM 




logM 


Z/Zo 




(Mo/yr) 


(10-*^*) 


(Mo/yr) 


2 


7.817 


-5.120 


8.511 


-5.815 


1/2 


8.554 


-5.631 


9.314 


-6.328 


1/5 


8.702 


-5.971 


9.474 


-6.665 


1/10 


8.751 


-6.225 


9.528 


-6.921 



well as the total amount of ionizing photons in the hydrogen 
and helium continua, agree within a few percent. The changes 
in the optical He i and He ii lines are also negligible. A decrease 
of at most 3% in equivalent width is observed compared to the 
constant scaleheight model. 

A fair comparison for the near-IR lines is more difficult, 
as many of these lines are partly formed in the transition zone 
where the photosphere connects to the stellar wind. In both the 
ISA-WIND as in the constant scaleheight model this connection is 
made in an ad hoc fashion. Keeping this in mind, we find that 
the differences in equivalent width between the two models is 
at most 4% for the He i Unes and at most 12% for the C iv Unes. 



2.2. Comparison with otiier model atmospliere codes 

Before comparing ionizing fluxes predicted by our models with 
those from other codes, we first briefly discuss the basic physics 
treated in the model atmospheres used for this comparison. 

We focus on a comparison of main sequence (dwarf) O star 
models comput ed using the cmfgen co de (this study), the WM- 
BASic code of Pauldrach et al.' and the CoStar models 

of Schaerer & de Koter (1997y. The latter make use of the isa- 
WIND code of Ide Koter et all JlOol Il997l) to predict the stel- 
lar spectrum. A schematic comparison of assumptions made in 
these three codes is given in Table 0] Note that a comparison 
of codes does not address the question concerning the level of 
realism of the spectra produced. Assessing this issue requires 
empirical testing by analyzing Hii regions containing single 
hot stars of which ideally the basic parameters have been deter- 
mined from detailed quantitative spectroscopy using the same 
models as those applied to pred ict the ionizing fluxes . This 
approach has been followed by ICrowther et alJ (Il999l) . who 
tested the cmfgen and isa-wind models for a nebula containing 
a cool WN star and found similar results except for the ion- 
izing flux distributions below the He i edge at 504 A, which 
were quite different. The next best thing is to study Hii re- 
gions containing single sta rs with we ll defined spectral types 
and effective temperatures. ISchaerej reviews the suc- 

cess of non-LTE model atmospheres in reproducing diagnostic 
line ratios of nebula irradiated by O-type stars. The metallicity 
dependence of the ionizing flu x may be studied by modelling 
extra-galactic Hii regions (e.g. 'Br esolin et al1ll999t loev et alJ 
HqOO; Bresolin & Kennicutt 2002). 

A comparison of ionizing fluxes from different codes is 
straightforward, however, interpreting the differences is - un- 
fortunately - not. Modest discrepancies are expected to origi- 
nate from small differences in the adopted abundances and pho- 
tosphere & wind parameters of the models one is comparing. 
For instance, a difference in effective temperature of 1000 K 
yields a flux difference in the Wien part of the spectrum of 10- 
15%. This may seem a needless 'mistake' and one that is easily 
fixed. However, these complex models are difficult to compute 
and it turns out to be less straightforward than one may think 
to tune all parameters. Having said this, the main sources of 
differences between codes hkely is in the way in which essen- 
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Table 4. Schematic comparison of assumptions made in the codes used for the EUV comparion. CMF implies a co-moving frame 
treatment of line transfer; SE denotes that the state of the gas is derived from statistical equilibrium (i.e. non-LTE). See text for a 
discussion. 



Treatment 


CMFGEN 


WM-BASIC 


CoStar 


Temperature structure 


Consistent 


Consistent 


Grey-LTE with blanketing feedback 


Density structure 








Photosphere 


Constant H 


Consistent 


Grey-LTE -I- p^-cont 


Wind 


j8-law 


Consistent 


yS-law 


Non-LTE treatment 


Consistent 


Consistent 


Consistent (except Fe) 


Line Transfer 


CMF 


CMF 


Sobolev 


Iron-group 


SE 


SE 


Modified Nebular 


Number of transitions/lines 








In SE 


-9000 


-30000 


-1000 


In spectrum 


-17 000 


-4000 000 


-100000 


Blanketing 


Consistent 


Consistent 


Blanketing factors 



tial physical processes are treated. Notably i) the treatment of 
the equation of motion, ii) the exact treatment of the non-LTE 
rate equations. Hi) the number of lines included (affecting the 
amount of line blocking), and iv) the way line blanketing ef- 
fects are taken into account. 

In the CMFGEN models we present here, we have described 
the photospheric density structure in a relatively simplistic way 
(see above). The CoStar models, using an atmospheric struc- 
ture calculated with the isa-wind code, have a more sophisti- 
cated density stratification. Here a grey LTE temperature struc- 
ture and an Euler equation taking into account continuum ra- 
diation pressure are iterated until convergence, accounting for 
the feedback of metal line blanketing on the Rosseland optical 
depth (see below). The wind is described in a similar way as 
for the CMFGEN models. The WM-basic models have the most 
self-consistent description of the density structure, as a final it- 
eration is done, based on the non-LTE temperature structure, 
in which both the continuum and line force is included in the 
equation of motion. This yields the density structure through- 
out the photosphere and wind, i.e. no ad hoc wind velocity law 
is required. 

The CMFGEN and WM-basic models compute the non-LTE 
state of the gas in a self-consistent way, using a super-level for- 
malism. In total these codes treat ~ 9 000 and ~ 30 000 transi- 
tions, respectively. In this respect, the CoStar models are rela- 
tively simple as the statistical equilibrium equations are solved 
explicitly for H, He, CNO and Si only, yielding - 1 000 tran- 
sitions, while the iron group elements, adding ~ 100000 lines, 
are treated in a modified nebular approximation. Moreover, the 
Sobolev approximation is applied to describe line transfer. Note 
that CMFGEN treats - 17 000 individual lines in the emergent 
spectrum, while WM-basic accounts for more than 4 000 000 
lines in the line-force and blocking calculations. 

Line blanketing - i.e. redistribution of flux due to line opac- 
ity and the feedback of this opacity on the temperature struc- 
ture as a result of the backwarming effect - is self-consistently 
treated in cmfgen and WM-basic, but is only accounted for in an 
approximate way in the CoStar code. To describe this in some 



detail: in CoStar, a Monte-Carlo formalism describes the line 
transfer problem of the metals. On the one hand, this allows to 
describe the feedback of blocking on the atmospheric structure 
using wavelength averaged blocking factors. On the other hand, 
this neglects photon conversion. This 'branching of photons' is 
relevant for redistributing flux, as it tends to 'soften' the (EUV) 
radiation field. 

The predicted EUV fluxes (especially in the He i contin- 
uum) are sensitive to line blanketing. In this respect, the cmf- 
gen and WM-BASic models include a more consistent treatment 
and may in principle be expected to yield more reliable results. 
For the ionizing flux, the prescription of the wind structure is 
relatively unimportant for main sequence stars, as these have 
modest winds. The treatment of blanketing and the specifica- 
tion of the photospheric density structure affect the strengths 
of lines, and therefore may have an effect on He i and He ii. 
As these lines are used to define spectral type, this may cause 
a modest uncertainty (up to one sub-type) in the spectral types 
assigned to the models. In this respect the WM-basic models 
are expected to behave best. 

3. The spectral energy distribution 

3.1. Comparison of predicted EUV fluxes 

Figure [2 shows a comparison between the spectral energy dis- 
tribution of comparable models computed with the three codes 
discussed above. The top left panel shows the emergent ion- 
izing spectrum for our CMFGEN model for reff=48 670 K 
(model #2). In the middle left panel, this cmfgen model is 
shown together wifli W M-BASic model D-50 (T^eff = 50000 K, 
'Pauldrach et al.''200l'^ and CoStar model C2 (T^ff = 48 529 K, 
.Scbaerer & de Koter 1997.) . The bottom left panel shows par- 
tially integrated q values of these three models. The spectra 
shown in the middel left panel are binned in wavelength to bet- 
ter bring out the overall flux behaviour, which is why they ap- 
pear smooth. However, a forest of metal lines does contribute 
to all of them. The Lyman jump is modest due to the high de- 
gree of ionization. From 912 A down to the He i jump at 504 A 
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Fig. 1. (Top panel) Emergent SED for models #2 (reff = 48 670 K) and #7 (T^g = 35 900 K). The continuum is plotted by a 
dashed line. (Middle panel) Comparison of the emergent SED for models #2 (left panel) and #7 (right panel) with CoStar and 
WM-BASIC dwarf models at similar T^g. (Bottom panel) The number of photons in cm"^ s ' below wavelength Aq calculated from 
the models considered in the middle panel is plotted as a function of this wavelength Aq. 
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(which is essentially absent), the most prominent features are 
strong resonance lines of N iv A765, O iv /1789, A609, A554 and 
O V A629. In this interval, the WM-basic model shows an on av- 
erage about 15% higher flux level compared to our cmfgen cal- 
culation. This can be explained by the temperature difference of 
the models. Taking Tetf effects into account, the CoStar model 
produces a somewhat softer spectrum. 

In the He i continuum the spectrum is dominated by a dense 
forest of mostly Fe iv-vi lines. Here the cmfgen and WM- 
BASic model show a roughly similar behaviour, both in terms 
of continuum flux and amount of line blocking. The CoStar 
model now shows a harder spectrum. The explanation for the 
CoStar predictions likely involves the treatment of the block- 
ing factors, which result from a Monte-Carlo simulation in 
which the effects of photon conversion (or 'branching') are 
neglected. This forces all radiative excitation/deexcitation pro- 
cesses to be scatterings, which prevents flux to be distributed 
towards longer wavelengths and therefore keeps the radiation 
field 'hard'. This would also explain the relatively soft flux in 
the CoStar spectrum at A > 500 A. Other effects that may 
play a (more modest) role are the temperature gradient in the 
He I continuum. A steeper gradient causes a higher flux at 
relatively short wavelengths (where the continuum is formed 
at relatively high Rosseland optical depth), while at relatively 
long wavelengths the flux will be lower. Also the neglect of 
line-broadening in CoStar models may have an effect (see 
ISchaerer & de Koter '1997"). Note that the exact location of the 
line forest differs between the cmfgen and WM-basic models. 
(To see this properly one first has to correct the continuum lev- 
els for differences in T^ff). This is most likely connected to dif- 
ferences in the predicted ionization of iron. The relevant iron 
ions each show a preferred region in which their line forest is 
concentrated. It holds very roughly that this is at 510 - 440 A 
for Fe iv; 460 - 320 A for Fe v; 340 - 260 A for Fe vi, and 
at 270 - 160 for Fe vii. This suggests that the dominant photo- 
spheric ionization of the WM-basic model is Fe vi, while Fe v 
is most important in cmfgen (see the pronounced difference at 
~ 400 A). 

Though not visible in the figure, the He ii continuum fluxes 
of the CMFGEN and WM-basic models agree reasonably well 
(which may be slightly surprising in view of the difference 
in T^etf)- Again the CoStar model produces a harder spectrum. 
The same reasons as discussed above may apply. However, it 
may also be connected to the use of the Sobolev approxima- 
tion and/or the velocity structure in the wind. The latter may 
cause a depopulation of the He ii ground state as a result of the 
velocity gradient in the wind, in combination with the temper- 
ature in the continuum forming layers at the wavelength of the 
He II resonance lines (Gabler et al. 1989.) . 

The right panels show a comparison between our cmfgen 
model for T^ff = 35 900 K (model #7), WM-basic model D-35 
having Tetf - 35 000 K, and CoStar model B2 for which the 
temperature is 36 300 K. Taking into account the differences 
in effective temperature, the EUV flux distributions compare 
well though again the Lyman flux at A > 700 A of the CoStar 
model is somewhat softer. Below ~ 700 A - i.e. again in the 
regime of severe flux blocking - the opposite behaviour is seen: 
the CoStar spectrum is harder. The likely reasons for this be- 



haviour have been discussed above. A detailed inspection of the 
energy distributions shows that the WM-basic models produce 
a clear drop-off in flux at the N iii and C iii ionization edges near 
260 A. The presence of this sizable bound-free jump is some- 
what surprising in view of the high T^^fi of the model, and is not 
seen in cmfgen and CoStar models. 

The fluxes in the He ii continua of the three models differ by 
orders of magnitude. Possible reasons for these discrepancies 
have been discussed above. 

We conclude that down to ~ 500 A the EUV flux distri- 
butions compare reasonably well, though the flux levels do 
differ. However, at shorter wavelengths large (order of mag- 
nitude) discrepancies may occur For the CoStar predictions, 
in comparison to the other two codes, the major part of these 
differences can be traced back to the treatment of blanketing. 
The models also show large differences for the Hen con- 
tinuum fluxes. This is a relevant problem with respect to the 
internal consistency of different code predictions. However, 
in O-type stars the dominant flux component at A < 228 A 
tends to be that of non-thermal soft X-ray emission arising 
from shocks in the s tellar wind (e.g. iMacFarlane etalJll994t 
iFeldmeier et al J1998). As no self-consistent description of this 
process is presently available, this seems a more acute problem. 

3.2. Comparison of ionizing fluxes 

Figure|2shows the total number of Lyman continuum photons 
emitted per second, Qo, and the total number per second per 
cm-, ^0 - QqI^ttRI, for both our O star main sequence models 
and the comparison models discussed above. Also given is the 
ratio between Lyman and He i continuum photons, 2o/2i (bot- 
tom panel). We added predictions from Kurucz models, pre- 
sented by Vacca et al. (1996); however, we will not include 
these in the discussion of the results. 

To eliminate differences in the adopted stellar radii, we first 
focus on the predictions for qo- The results are plotted versus 
Teff, allowing us to assess the differences in as a result of 
small differences in adopted temperatures. All predictions are 
very similar, only showing discrepancies up to 0. 1 dex. Why is 
this so? 

The principal answer is that though line blanketing effects 
are strong at wavelengths shorter than ~ 500 (700) A in the case 
of model #2 (#7), the flux is mainly redistributed within the 
Lyman continuum. Only limited energy 'leaks' out at A> 912 
A. This may be illustrated using a simple estimate. The line for- 
est removes ~ 40% of the flux in both model examples. These 
photons are back-scattered and re-thermalized. Assuming the 
thermal emission may be characterized by the emerging con- 
tinuum flux distribution, one expects 60% (20%) of this en- 
ergy to be emitted still within the Lyman continuum. Therefore, 
only 16 % (32%) leaks out into the Balmer continuum and be- 
yond. In terms of number of ionizing photons, the 'damage' 
is even less as redistribution of energy towards longer wave- 
lengths within the Lyman continuum requires an increase in 
the number of photons. This explains why the number of ioniz- 
ing photons only differ up to ~ 0.05 and 0. 1 dex for models #2 
and #7 respectively. The neglect of branching in CoStar mod- 
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Fig. 3. (Top panels) Comparison of the EUV spectrum of solar abundance model #5 (dotted line) with equivalent models with 
Z - I / 10 Zq (left panel) and Z - 2 Zq (right panel), plotted in grey. (Bottom panels) Binned version of the spectra plotted in the 
top panel. 



els is separate from the above discussion (it concerns forward- 
scattered photons). However, note that the iron forest is dom- 
inated by transitions from meta-stable levels. The level struc- 
ture of the relevant iron ions is such that branching typically 
produces two photons, one at A < 912 A and one at A > 912 A. 
Therefore, also this effect preserves Lyman continuum photons, 
though it does affect the hardness ratio Qi/Qq. The former is 
illustrated in the bottom panels of Fig. [2 where we show par- 
tially integrated q values. 

The hardness ratio Qq/Qi, also independent of stellar ra- 
dius, yields agreement within ~ 0.2 dex. The CoStar models 
having Teff > 40000 K show a harder spectrum than do the 
CMFGEN and WM-BASic models - as discussed above. This is in 
agreement with the flux differences in the Lyman and He i con- 
tinua discussed for model #2. 

The apparent reasonable agreement between number of 
ionizing photons is somewhat misleading. Nebular lines from 
ions such as Ne iii or S iv require radiation from below 500 



A. In this spectral regime, as we have seen above, the ionizing 
fluxes may differ strongly implying large differences in emis- 
sion line strength. 

3.3. EUV fluxes and Qq and Qi values as a function of 
metal content 

In Table |5] we show the predicted number of ionizing photons 
Qo and Q[ for models #2 and #5 as a function of metal abun- 
dance. The results show that a decrease in Z of a factor of 
twenty only alters the number of Lyman continuum photons 
by 0.03 dex. Over the same metallicity range the number of 
He I continuum photons increases by up to 0.08 dex. 

The modest dependence on metal content is again the result 
of the way line blanketing redistributes the flux (see above). In 
the top left panel of Fig.|3lwe show a comparison of the EUV 
spectrum of solar abundance model #5 with that of a model 
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Fig. 2. Variation of the Lyman continuum ionizing photons Qq 
(top panel), qo (middle panel), and the hardness ratio Q\/Qo 
(bottom panel) with effective temperature. The results from 
our CMFGEN models (solid line) are compared to those from the 
WM-BASic models (dotted line), CoStar models (dashed line) 
and Vacca et al. ( 1996) (dash dotted). We use WM-basic mod- 
elsD-30, D-35, D-40, D-45, D-50 and D-55 ( Pau ldrach et all 
[2001'), and CoStar models A2, B2, C2, D2, E2 and F2 
(Schaerer & de Koter 1997). 

Table 5. Predicted number of ionizing photons go and Q[ for 
# 2 and #5 as a function of metal abundance. 



Model #2 



Model #5 



Z/Zq 


loggo 


loggi 


logQo 


loggi 


2 


49.70 


49.16 


49.09 


48.40 


1 


49.69 


49.18 


49.08 


48.44 


1/2 


49.69 


49.19 


49.07 


48.46 


1/5 


49.68 


49.20 


49.07 


48.47 


1/10 


49.67 


49.21 


49.06 


48.48 



with identical basic parameters, only with Z = I/IOZ©; in the 
top right panel the solar and twice solar model are compared. 
Relative to the 1/10 Zq model, the solar metallicity prediction 
clearly shows an increased flux level in the region from the 
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Fig. 4. The effect of the metal content on the number of pho- 
tons in cm"^ s"' below a certain wavelength Aq is plotted as a 
function of this wavelength Aq. The effect is shown for models 
#2 and #5. 

Lyman jump down to ~ 600 A, which is compensated by a 
decreased flux below 600 A. This is best seen in the bottom 
left panel, where we show binned versions of the spectra. The 
tendency to redistribute flux within the Lyman continuum is 
also illustrated in Fig.|4] where we show partially integrated q 
values. 

The conservation of ionizing photons works less well for 
Qi, as expected. Here a decrease in metallicity of a factor 
twenty results in up to 20% more ionizing photons. This im- 
plies that those nebular species that depend on a specific EUV 
flux range below ~ 500 A may show a strong dependence on 
metallicity. We will discuss this in Sect.ls] 

The results for the number of hydrogen ionizi ng photons 
are similar to those obtained by Smith et al. ('2002') using uni- 
fied WM-BASic models, by Lanz & Hubenv ( 2003) using plane 
parallel TLUSTY models and bv iKudritzkil J2Q0A . who used 
WM-BASic models of very massive O stars at Z = 10""* Zq to 
Z = Zq. The results of the two former studies show that the 
convervation of Qo for dwarf stars breaks down at lower effec- 
tive temperatures (at 32 000 to 37 500 K depending on model 
assumptions). Compared to our mo dels the ri umber of He i con- 
tinuum photons in the models of Smith et alJdo not show a si m- 
ilar behaviour The models of Kudr itzkil and lLanz & Hubenvl do 
show a conservation of the number of photons, where the lat- 
ter models again show that this convervation breaks down for 
lower effective temperatures (< 45 000 K). 

4. Spectral classification 

4.1. Optical classification 

The optical spectral classification of O stars relies on a 
quantitative criterion based on the ratio of the equivalent 
width s of He: /1447 1 and Hen ^4542 fconti & AlschuleJ 
.197 It IConti & Fros^ IT977t Maflivs. .1988.) . In 09 stars, flie 
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He I /14471 line is stronger; the ratio changes gradually 
so that the Hen /14542 line begins to dominate at type 
06 and earlier. According to this criterion, the value of 
log[W^(4471)/W^(4542)] unambiguously determines the spec- 
tral type of the star. 

This observational criterion can be used to assign spectral 
types to our models. Fig.lsjcompares the equivalent widths of 
the He I and Hen classification lines observed in dwarf stars 
with the predictions from our models. The slopes of the dot- 
ted lines correspond to the ratios of the equivalent widths de- 
limiting the O star subclasses according to the calibration of 
[Mathvs ( 198^. We investigate the assigned spectral classifica- 
tion for the eff'ects of varying respectively log g, metallicity and 
microturbulent velocity. Before discussing these dependencies 
in detail, we will first elaborate on how well our initial model 
grid for dwarf stars fits these observations. 

The observational da ta for the dwar f stars are taken from 
IConti & AlschuQl[l9Ti])jConti.(.197i.lConti & Frost (1977l 
and iMathvsl lll988l Il989l) . The lumir iosity class wa s veri- 
fied using the spectral classification of 'Walborn| (ll97lL fl972l 
Il97d. .1973. .198Z) . Stars that were not classified as dwarfs by 
Walborn were discarded. S pectroscopic binaries w ere identi- 
fied using the catalogue of Pedoussaut et al. ( 1996) and stars 
with variable radial velocity using Mathvs (.1988, 1989) . These 

were also removed . Finally, variable stars were deselected us- 
j ' ' 1 1 I 

ing the catalogue of lKholopov et al. ( 1998). In the top left panel 
of Fig.|5lwe show the selected observed data (using open di- 
amond symbols), including the error bars. The observations 
cover a strip in the Hei AAAll versus Hen AA5A2 equivalent 
width diagram that is broader than can be accounted for by un- 
certainties in individual measurements. The width of this strip 
therefore also reflects differences in stellar properties such as 
gravity, mass loss, metallicity and/or photospheric properties 
such as microturbulent velocity. 

4.1 .1 . Mass loss 

Also plotted in the top left panel are the models listed in 
Table ^ The predictions reproduce the observed range of 
He I AAAl 1 equivalent width values well. The equivalent width 
of He II AA5A2 seems to be overestimated by approximately 10 
to 20% for early- and mid-O spectral types. A possible expla- 
nation for this overprediction may be an adopted mass loss that 
is systematically too low. A higher mass loss would reduce the 
equivalent width, as wind emission fills in m ore of the photo- 
spheric absorption profile (Lanz et al. 1996). This effect is es- 
pecially important for the earliest spectral types, as these have 
the strongest winds. The top right panel of Fig. |5l shows the 
effect of an increase in M. Here a grid with doubled mass loss 
rate, indicated by squares, is compared to the standard grid (cir- 
cles). A change of this magnitude is charact eristic of the typ- 
ical error in both mass-loss determinations jPuls et al.lll996l) 
and mass-loss predictions ( Vink et al. 20001) . The line strength 
of He II AA5A2 is decreased by up to 10 to 20% in the early type 
models, as a result of wind emission. A modest decrease in the 
He I AAAll strength is also observed. 



4.1.2. Gravity 

The width of the observed strip of equivalent width ratios is 
partly explained by gravity effects. Kudritzki et al. ( 1983) al- 
ready showed that the ratio W.^(4471)/W^(4542) depends on 
luminosity class. Here we show that gravity variations within a 
luminosity class also significantly effect this ratio. The main se- 
quence gravities adopted by Vacca et al. (1996) should be seen 
as mean values for the ensemble of stars of given spectral type 
used for their calibration. We have calculated a grid of models 
in which we increased and decreased the stellar gravity by 0.3 
dex. This roughly corresponds to the typical difference between 
dwarf and giant stars. The bottom left panel of Fig.|5]shows the 
predicted strip of equivalent width ratios spanned by this grav- 
ity range. Note that e.g. a decrease (increase) of log g by 0.3 dex 
implies about one spectral sub-type earlier (later). A decrease 
of log g implies a decrease in electron density, which causes 
a decrease in the recombination rate, effectively increasing the 
helium ionization. The opposite is true in case of an increased 

gravity. 

The new temperature scale of O stars (iMartins et al."2QQ^ 
see also de Koter et al. 1998) implies a decrease in Teff of about 
10 percent compared to the calibration of IVacca et alJ (Il996l) 
that is used for our grid. The lower luminosity and mass im- 
plied by this new Teff results, however, in only a modest shift of 
~ 0.06 dex in gravity. Therefore, this temperature recalibration 
can only account for a very minor part of the overprediction of 
He II ^4542. 

In part the overprediction of He ii AA5A2 may also be the re- 
sult of a to high stellar mass. Indeed, a systematic difference has 
been noted in masses that are based on spectroscopic analysis 
of Balmer line profiles an those that result from a comparison 
of basic parameters (L, and T^s) with evolutionary tracks (e.g. 
Herrero et al. 1992). Recent new results (Herrero et al. 2Q02) 
appear to show that most of this discrepancy can be explained 
by improved spectroscopic modeling and the lower T^^ val- 
ues for O-stars.^illj adopting the spectroscopic masses as pre- 
sented bv lVacca et al] (fl996l) (which are 30 to 70% less for the 
earliest spectral types) yields a decrease in He ii AA5A-2 equiv- 
alent width of about 10 to 15% for the earliest spectral types. 
This should be seen as an upper limit of a possible mass effect. 

4.1 .3. metallicity 

metallicity effects may also explain part of the scatter in equiv- 
alent widths, but are not expected to be responsible for the 
systematic shift in the Hen AA5A2 line strength. Variations 
in the stellar metallicity are undoubtedly present in Galactic 
O stars. Emission-line studies of H ii regions, planetary neb- 
ulae and supernova remnants consistently show that the aver- 
age abundance of oxygen decreases with Galactocentric dis- 
tance by ~ 0.06 dex per kpc, with values of ~ 2.5 Zq in the 
Galactic Center, and down to ~ 1/4 at Galactocentric dis - 
tances larger than 15 kpc (see e.g. iHenrv & Worthev"l999t). 
Low abundance systems such as the Small Magellanic Cloud 
have metallicities of the order of ~ 1/8 Z0. The top right 
panel in Fig. |5] shows the variation of the equivalent widths 
of He I /1447 1 and He n /14542 as a result of changing the atmo- 
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Fig. 5. Top left panel compares the equivalent width of Hei A4471 and Hen /i4542 measured in the model grid with observed 
values in main sequence (dwarf) stars (open diamonds). Bottom left panel shows the effects of an increase (triangles) and a 
decrease (squares) of 0.3 dex in log^. In the top right panel a grid with doubled mass-loss (squares) is compared to the model 
grid (circles). Also shown are the equivalent widths of models #2 and #5 with Z - 2, 1/2, 1/5 and I/IOZ0 (triangles). The 
decrease in metallicity results in an increase of the equivalent width of He 1 4471 A. The bottom right panel compares a grid with 
Vtuib = 5 kms"' (circles) to a grid with Vmrb = 20 kms ' (triangles). The dotted lines indicate the parameter space delimited by 
every subclass according to the empirical classification bv .Mathvs.(.1988.) . 



spheric metallicity - and its indirect effect on line blanketing 
and stellar wind mass loss - from 1/10 to 2 Zq. The models 
with modified metal content are indicated by solid triangles. 

Two effects related to metal content may influence the line 
strength of the Hei AAAll and Hen AA5A2. First, an increased 
metallicity implies an increased line blanketing and therefore 
stronger diffuse field. This will increase the ionization of he- 
lium. The top right panel of Fig.|5]indeed shows this behaviour 
for model #5, i.e. a decrease in the He i and an increase in the 
He 11 line strength. Second, an increased metallicity implies a 
stronger stellar wind, therefore an increased filling in of line 
profiles by wind emission (see also e.g. Herrero.2003i) . This ef- 



fect is more important for earlier type stars, which have a larger 
mass loss. This second effect dominates in model #2, resulting 
in a decreased He 11 AA5A2 line strength with increased metal- 
licity and vice versa. 



We conclude that metallicity effects may change the spec- 
tral typ e by up to one subtype. This is in agreement with find- 
ings of lPuls et al.l (l2003). which showed that the inclusion of 
line blanketing at SMC metallicities resulted in a much more 
modest change in Teft relative to that found for Galactic stars. 
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Fig. 6. Effective temperature shift with spectral type between 



the lVacca e t al. ( 1996) calibration and the temperature scale de- 
rived from the models with Z - Zq and Vmrb = 5 km s ' . 



4.1.4. Microturbulence 

The effect of microturbulence is expected to explain part of the 
scatter in the observed He: /14471 and Hen A454-2 equivalent 
width ratios and is shown in the bottom right panel of Fig|5] The 
standard grid assumes a constant microturbulent velocity of 
5 km s"' . We calculated a new grid adopting Vtuit = 20 km s"' , 
producing stronger lines. A comparison shows that the Hen 
line is almost insensitive to the introduced variation in micro- 
turbulent velocity. The saturated He i line however shows a sig- 
nificant change in strength for models of spectral type 05 and 
later. The relative changes are 4% for the model with the weak- 
est He I line and 12% to 20% for the models with the strongest 
He I lines. An increased microturbulence may shift the assigned 
spectral type by up to half a sub-type towards later spectral 
type. The impact of microturbulence in the line formation cal- 
culations of H and He i in OB stars, and consequently, in the de- 
rived stellar atmosph eric parameters, has been investigated by 
several authors (e.g. 'Smith & H owarthlll998t iMcErlean et alJ 
[l99 8: Villamariz & Herrero 200^! 

Summing up the above discussion, we conclude that all-in- 
all the predicted equivalent width of He i A44-1 1 and He ii A4542 
reproduce the observed behaviour reasonably well, though the 
Hen line strength is somewhat overpredicted for early- and 
mid-O types. We also show that the spectral type of dwarf stars 
is not only a function of effective temperature of the star, but 
also depends, albeit to a lower degree, on microturbulent veloc- 
ity, metal content and on gravity, i.e. gravity variations within 
the dwarf class. Variations in these parameters explain the ob- 
served strip of equivalent width of He i A441 1 and He ii A4542 
for given spectral type. 

Finally, when comparing the spectr al types of our mode ls to 
the original spectral types assigned bv lVacca et alJ ( Il996l) . we 
see that there is a systematic shift to earlier spectral types for 
given temperature. In other words, the T^ff scale derived with 
our models shows a systematic shift to lower temperatures. 
Fig-Elshows the difference for Z = Zq and Vtuit = 5 kms"'. 
The decrease from the iVacca et alJ scale to our new scale varies 



between ~ 2500 K and ~ 3800 K and is largest for stars in the 
spectral range of 05 to 09. This decrease of T^ff due to line 
blan keting effects h as recently been found in several studies 
re.g. Ide Koter et al.lil 998t iCrowther et al.i.2o'oil iMartins et al J 



4.2. K-band classification 

K-band (2 fj.m) spectroscopy has been used in the past few 
years to identify and classify the characteristics of newly- 
form ed O stars that are found in heavily obscured H ii regions 
(e.g. IWatson et alJll997l IWatson & Hansonlll997i iKaner et alJ 



'2002'. 'Bik et al.'2003'. "M artin-Hernandez et al.l2003l) . The typ- 
ical dust extinction of more than 10 mag at visual wavelengths 
towards these Hn regions has restrained optical photometry 
and spectroscopy of such stars. At near infrared wavelengths, 
the extinction is much less, and direct observations of their pho- 
tospheres are possible. At longer wavelengths, observations fail 
to detect their photosphere directly because the spectral energy 
distribution is dominated by emission from the surrounding 
dust. Hence, it becomes essential to understand the behaviour 
of the photospheric features in the near-IR regime, the only ob- 
servational window of such young O stars. In this section we 
compare the K-band spectral behaviour of our models with ob- 
servations. 

iHanson et alJ ( 1 19961) present an atlas of K-band spectra of 
a large number of relatively nearby, well-studied, optically vis- 
ible massive stars with the intention of investigating the varia- 
tion of 2 fj.m spectral features with spectral type and luminos- 
ity in known OB stars. They find that in early-O dwarf stars 
(03, 04, some 05), Hen /12.189 (wavelengths of the near- 
IR lines are in micrometers) absorption and N iii (observed at 
2. 116 jum) emission are the dominant spectral indicators. Mid- 
O stars (05, 06, some 07) show C iv /12.078 emission and can 
begin to show He: /12.113 absorption. Late-O (08, 09) and 
early-B stars lack He 11, N iii and C iv, showing instead He i and 
fairly strong (in comparison to hotter O stars) Hi /12.1661 and 
Brackett y. 

Our grid of models can be used as a first test against the 
features found in the K-band spectrum of O-type dwarf stars. 
Many complications are to be expected in the modeling of 
these lines since they are formed in the transition region from 
the stellar photosphere to the super-sonic stellar wind (see e.g. 
iKudritzki APuls"2000^. A full parameter study is highly desir- 
able to understand their diagnostic value. In this paper, we only 
discuss the temperature and microturbulent velocity effects and 
refer to Lenorzer et al. (in prep.) for a detailed parameter study. 

Figure shows, as an example, the comparison between 
the synthetic K-band spectrum from model #4 with Z = Zq and 
Vturb = 10 km s"' (which corresponds to an 05 main sequence 
sta r) and the spectrum o f HD 93204, an 05 V((f)) star observed 
bv lHanson et alJ lll996l) . The synthetic spectrum was degraded 
to the resolution of the observed spectrum (A/AA ~ 1500). 
The He I lines at 2.058 and 2.113 fj.m, Bry, and the Hen 
line at 2.189 fj.m are clearly visible in the synthetic spectrum. 
Qualitatively, the behaviour of these lines closely resembles 
that shown in the 2 //m atlas. He i /12.058 is almost undetectable 
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Fig. 7. The synthetic K-band spectrum from model #4 with 
Z = Zq and Vtmb = 10 kms"', which corresponds to an 05 
main sequence star (dashed line) is compared to the spectrum 
of the 05V((f)) star HD 93204 (solid line). The synthetic spec- 
trum is degraded to the resolution of the observed spectrum 
(/i/A/l ~ 1500). The lines observed in the K-band are indicated. 
Note that while the C iv doublet transition is clearly observed 
in emission, our model predicts it to be in absorption. The N iii 
emission at 2.1 16 /im is not included in the models. 



in our three hottest models but starts to show a strong absorp- 
tion in the interm ediate tem perature models to finally weaken 
in model #8. .Hanson et alJ tl996) observed that this line may 
turn into emission for B stars. Bry and Hei /12.113 strengthen 
towards later-type stars, while Hen /i2.189 becomes progres- 
sively weaker 

The models show a discrepant behaviour in the case of 
the metal lines. The Nm emission at 2.116 jjm, with transi- 
tion 8-7, is simply not included in the models. The C iv dou- 
blet t ransition (3pP°-3d^p ), observed in emission 04 to 08.5 
stars ( iHanson et al.ll 19961) . is accounted for in our models. It 
is, however, predicted in absorption in all models that have 
Te« ?J 41 000 K, i.e. spectral type 05 or earlier. The disap- 
pearance of these lines in late O-type stars is well predicted 
and is the result of the recombination of C iv to C iii. The ori- 
gin of the observed emission lines is more difficult to under- 
stand. Recent high spectral resolution (R ~ 8000) observations 
of O-type stars (Bik et al. in prep.) resolve these C iv lines in 
a few stars, yielding a FWHM of about 40 kms ~^ Moreover, 
the strength of the lines (e.g. Bry) as observed bv lHanson et alJ 
is not found to be much stronger in supergiants than in dwarf 
stars. These arguments strongly suggest that the near-infrared 
C IV lines have a predominantly photospheric origin, therefore 
their emission cannot be attributed to an enhanced C iv abun- 
dance and/or wind effects. If the emission was due to a temper- 
ature inversion in the line forming region, this would similarily 
affect He i /12.058. Indeed, both lines have regions of formation 
that largely overlap such that a temperature inversion would 
also be visible in the He i line. Our models, however, reproduce 
the He i line reasonably well. The transition of C iv observed 
in the near-infrared spectrum of O-stars involves low atomic 
levels that are populated via transitions located in the UV and 
far-UV. It is therefore likely that we do not reproduce the be- 
haviour of these lines due to a mistreatment of this complex 
spectral region. This issue is still being investigated but is en- 
couraged by recent calculations including C iii which show the 
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Fig. 8. Comparison between observed equivalent widths of 
Bry, He I /t2. 058 and Hen A2.1S 9 in main sequence stars 
(open squares, !Hanson et al ]ll996^) and those measured from 
our models with Vturb - 5 (solid line), 10 (dotted line) and 
20 kms '(dashed line). 



C IV doublet transition in emission. We therefore suspect that 
line pumping is the cause of these emission lines (Lenorzer et 
al. in prep.) 

We carried out a quantitative comparison with observed 
equivalent widths using the Hei /12.058 , Bry and He ii /12.189. 
Fig-IHcompares the predicted equivalent widths of these three 
lin es adopting Vfn^h = 5^ 10 and 20 km s ' , with those measured 
bv lHanson et alJ(ll99 6V Spectral types of the models are those 
derived from the optical lines (see Fig.|5}- The microturbulent 
velocity has no significant effect on the lines. The models fit 
reasonably well the observed equivalent widths, although they 
overestimate the strength of Hen /12.189 for spectral types 03 
and 04. 

We can investigate whether the equivalent widths of the 
2.058 fim Hei, Bry and 2.189 //m Hen lines are affected by 
variations in the metal content. For model #2 with Vmrb - 
10 kms"' (which corresponds to an 03 star), a change in 
metallicity from 2 to 1/10Zq produces a decrease in the equiv- 
alent width of the Hen line from 3.3 to 1.8 A. Bry shows only 
variations at the 30% level and the He i line is practically ab- 
sent at all metallicities (equivalent width < 0.1 A). In the case 
of model #5 with Vtmb = 10 kms ' (an 06 star), the equivalent 
width of the He i line remains practically unchanged; the Bry 
and He n fine show variations of the order of 20-30%. 
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Summing up, variations in microturbulent velocity have no 
significant impact on the equivalent width of the Bry, He i and 
He n lines, while variations in metal content have an apprecia- 
ble influence on the equivalent width of the He ii line in the case 
of very hot stars. 



5. The ionizing structure of single star H ii regions 

Observations of H ii regions combined with detailed photoion- 
ization models can be used to test the EUV spectrum of O 
stars. In this respect, the Infrared Space Observatory (ISO, 
iKessler et al. 1996) provided an unique opportunity by measur- 
ing the i nfrared (2.3-196 um) spectra of a large sample of H ii 
region s dOiveon et al.ll2002t IPeeters et af]l2002t IVermeii et all 
I2OO2I) . giving access to four elements (N, Ne, S and Ar) in 
two dififerent ionization stages. Ratios of fine-structure lines 
such as [Niii]/[Nn] 57/122 fim, [Ne iii]/[Ne n] 15.5/12.8 //m, 
[Siv]/[Sra] 10.5/18.7 jum and [Ariii]/[Arii] 9.0/7.0 //m probe 
the ionizing stellar spectrum between 27.6 eV (450 A) and 
41 eV (303 A) and can be used as indicators of the realism 
of stellar models at these energies. 

Basically, the ratio of two successive stages of ionization 
X^' and X^'^' of a given element X indicates the state of ion- 
ization of the nebula, which depends principally on the shape 
of the SED, more specifically on the number of photons able 
to ionize X"^' compared to that of Lyman continuum photons 
dVflchez & PageL1988 ). and on the ionization parameter U. In 
the case of an ionized sphere of constant gas density, n, and fill- 
ing fact or, e, U is defined by U = O n/(4nR^nc) oc (Qone^)^'^ 
(see e.g. Stasinska & Schaereiil997h . where is the radius of 
the Stromgren sphere and c is the speed of light. Any combi- 
nation of Qq, n and e that keeps the product Q^ne^ constant 
results in a similar ionization structure. 

In order to investigate the dependence of the fine-structure 
line ratios on the SED, nebular models have been computed 
with the photoionization code CLOUDY version 96.00-beta4' 
using MICE, the IDL interface for CLOUDY created by H. 
Spoon^. We computed nebular models for static, spherically 
symmetric and homogeneous gas distributions with one ioniz- 
ing star in the center. An inner cavity with a radius equal to 
10'^ cm is set, while the outer radius of the Hii region is de- 
fined by the position where the electron temperature reaches 
1 000 K. The grid of stellar models based on the cmfgen code 
and described in Sect.|2lis taken to describe the incident SED. 
We take models with Vtuit = 10 kms"'. The stellar spectra 
were rebinned to 2 000 points, which is the limit of points for 
the input SED accepted by CLOUDY. The nebular metallicity 
was set equal to that of the stellar model used to calculate the 
SED. Finally, with the aim to explore the dependence of the 
fine-structure line ratios on the SED alone, we decided to fix 
the ionization parameter U, so that the variations of the line ra- 
tios are only due to variations of the SED appearance. Hence, 
we fix Qq to the typical value of 10^** photons s"', the density 
to n = 10^ cm"^, and the filling factor to e = 1. We stress here 



' see http://thunder.pa.uky.edu/cloudy/ 

^ see .http://www.astro.rug.nl/~spoon/mice.htmII 



that unless U is constrained, these line ratios cannot be used to 
derive an absolute value for T^s (Stasinska & Schaerer 1993- 

We computed several sets of nebular models. First, a set at 
solar metallicity where the dependence on effective tempera- 
ture is analyzed, and second, a set where the effective temper- 
ature is fixed and the metal content is modified. This last grid 
allows us to investigate the dependence of the fine-structure 
line ratios on metallicity. 

Figure |9l presents the correlations observed between the 
fine-structure line ratios [Arm]/[Arii] 9.0/7.0 jum, [Siv]/[Sin] 
10.5/18.7 /im a nd [Neml/rNenl 15 . 5/12.8 ^m for the sam- 
ple of Galactic dOiveon et al .' 2002; Martm-Hernandez et al.' 
l2002a) and Magellanic Cloud H 11 regions ( V ermeii et al..2002.) 
observed by ISO. The ionization potentials of Arii, Siii and 
Ne II are, respectively, 27.6 eV (450 A), 34.8 eV (357 A) and 
41.0 eV (303 A). We do not include the ratio [Niii]/[Nii] 
57/122 yum in this figure; the ionization potential of Nii (29.6 
e V) is close to that of Ar 11, and hence, this ratio roughly probes 
the same range of energies as [Arm]/[Arn]. The data span 
a range in ionization of more than two orders of magnitude. 
Interestingly, the Magellanic Cloud H 11 regions nicely overlap 
with the trend o bserved for the Galactic objects at the high ion- 
ization end (see lMartm-Hernandez et alJ2002bl) . 

The diagrams in Fig. |9] probe the shape or hardness of 
the stellar spectrum between 303 and 357 A (panel a), 303 
and 450 A (panel b), and 357 and 450 A (panel c). In this 
diagram the results for the set of nebular models at solar 
metallicity are confronted to the observations. As a com- 
parison, we also show the results of grids of models com- 
puted using the WM-basic and CoStar stellar models. We use 
WM-BASic models D-3 0, D-35, D-40, D-45, D-50 and D-55 
(|Paiildrachetal.ll200lh. and CoStar m odels A2, B2, C2, D2, 
E2 and F2 ("Sc haerer & de Koteill 19971) . For all model sets, the 
predicted nebular excitation increases with increasing effective 
temperature of the ionizing star. However, differences among 
the models, and between the models and the observations, are 
clearly present. The largest differences are seen in panels (b) 
and (c), for which the range of energies tested is larger. Several 
conclusions can be gleaned from the inspection of Fig.|9] (1) 
For the same nebular conditions, i.e. same ionization param- 
eter U, the nebular models based on the CoStar SEDs with 
Teff > 36 000 K largely overpredict [Ne iii]/[Ne 11] for a given 
[S iv]/[S III] or [Ar iii]/[Ar 11] when compared to the observa- 
tions and the other SEDs. This can be understood in terms of 
an excess of photons below 303 A. This effect is evident in 
Fig.[0 which clearly shows that the CoStar stellar spectra are 
much harder at these energies than the other two atmosphere 
models. (2) The nebular models based on the WM-basic SEDs 
show a too soft stellar spectrum between 303 and 450 A for the 
models with T^s < 45 000 K, which is reflected in an under- 
prediction of the [Ne iii]/[Ne 11] and [S iv]/[S m] line ratios for 
a given [ Ar iii]/[ Ar 11] . This is also illustrated in Fig.^| which 
shows the variation of the [Ne in]/[Ne 11] and [Ar m]/[Ar 11] line 
ratios with T^s. The [ Ar iii]/[ Ar 11] line ratio, which is sensitive 
to the number of photons able to ionize Ar II (below 450 A) 
relative to the total number of Lyman photons, is, for a given 
U, practically independent of the stellar model used in the pho- 
toionization model. However, large discrepancies appear in the 
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Fig. 9. Comparison of the infrared fine-structure line ratios [ Ar iii]/[ Ar ii] 9.0/7.0, [Siv]/[Sra] 10.5/18.7 and 
[Ne iii]/[Ne ii] 15.5/12.8 //m observed in Galactic (open squares) and Magellanic Cloud (open triangles) Hn regions with 
predictions from photoionization models at solar metalUcity based on three different SEDS: cmfgen (solid circles), WM-basic 
(solid diamonds) and CoStar (solid triangles). A unique ionization parameter, U, is adopted for all the models. We stress that 
unless U is constrained, these line ratios cannot be used to derive an absolute value for T^g. 



case of the [Ne iii]/[Ne n] line ratio. This is because the range of 
energies probed by [Ne iii]/[Ne n] depends greatly on the way 
the photosphere and wind parameters of the stellar atmosphere 
are defined. Hence, the large differences between the models. 
In particular, the nebular models based on the WM-basic SEDs 
predict a lower [Ne m]/[Ne ii] than the other sets of models for 
Teff < 45 000 K. (3) In the case of the nebular models based on 
the CMFGEN SEDs, the discrepancies with the observations are 
smaller than for the other two model grids, thus implying that 
the stellar spectrum between 303 and 450 A is better described 

by CMFGEN. 

Finally, we quantify how variations in metallicity modify 
the ionization structure of Hii regions and hence, the inter- 
pretation of the above fine-structure line ratios in terms of the 
stellar content of the nebula. Observationally, the [Ne iii]/[Ne ii] 
line ratio observed i n H n regions show a loose corre lation with 
nebular metallicity jMartfn-Hernandez et alJ2002bl) . Figure fTTI 
shows the eff'ect of varying the metallicity for models #2 and 
#5 on the [Ariii]/[Arii] and [Ne iii]/[Ne ii] line strength ratios. 
In Sect. l3.3l it was discussed that, for a given eff'ective temper- 
ature, the stellar spectrum softens with increasing metallicity. 
This effect has a clear impact on the degree of ionization of the 



H n region: the fine-structure line ratios decrease with increas- 
ing metallicity. The biggest effect occurs for the [Ne iii]/[Ne ii] 
ratio, which gets reduced by about a factor of 10 when the 
metallicity increases from 1/10 to 2Zq. Hence, the effect of 
metallicity is such that, for instance, the ionization structure as 
traced by [Ne m]/[Ne ii] of a half solar metallicity nebula ion- 
ized by a star with T^^ - 41 010 K is, for the same ionization 
parameter, almost identical to that of a solar metallicity nebula 
with a central star of T^s - 48 670 K conform Fig.[^(also see 
iMorisset et al.l2003h . 



Infrared line ratios have been used to infer the stellar con- 



tent of H II regions (e 




'Takahashi et al." 


'2000^, 


"Okamoto et alJ 


.2001: 


Morisset et al. 


2002), and to constrain the ages in 


starbursts (e 




ICrowther et alJ Il999l iThornlev et alJ 1200(1 


ISooon et al."2000V However, the laree 


dependence of these 



fine-structure line ratios on metallicity make them unsuitable 
for such studies, unless the metallicity is known, and more- 
over, its influence on the stellar energy distribution is properly 
accounted for. 



M. R. Mokiem et al.: Metallicity and the spectral energy distribution and spectral types of dwarf 0-stars 



15 



0) 



W) -1.0 - /, 



M 
O 




• CFMGEN 

♦ WM-basic 

A CoStar 



35 



40 45 

Terr (10' K) 



50 



55 



Fig. 10. Predicted [Ariii]/[Arii] 9.0/7.0 //m and [Nera]/[Neii] 
15.5/12.8 i^m fine-structure line ratios as a function of T^ff for 
the three sets of models at solar metallicity considered in the 
text. A unique ionization parameter, U, is adopted for all the 
models. We stress that unless U is constrained, these line ratios 
cannot be used to derive an absolute value for Teff . 

6. Conclusions 

We have presented a grid of eight main-sequence (dwarf) star 
models ranging in effective temperature from 33 000 up to 
5 1 000 K. The models were constructed using the cmfgen code 
of lSillier & Miller. (J 998ii) . In order to investigate the variations 
of the spectral appearance on metallicity, and hence, their im- 
plications for spectral classification and the ionization structure 
of H II regions, we have calculated additional sets of models in 
which we varied the metal content from 2 to 1/1 OZ0. The main 
conclusions of this study are: 

• The total number of Lyman photons emitted is found to be 
almost independent of line blanketing effects and metallic- 
ity for a given effective temperature. This is because the 
flux that is blocked by the forest of metal lines at A < 600 
A is redistributed mainly within the Lyman continuum. 
Though some flux is removed from the ionizing continuum, 
the softening of the spectrum helps to conserve the number 
of ionizing photons. 

• We investigate the influence of metallicity, microturbulent 
velocity and gravity g on the optical lines used to spectral 
type O stars: He i A4411 and He ii /14542. We find that spec- 
tral type, as defined by the ratio of the equivalent widths of 
these lines, is not only a function of the effective tempera- 
ture of the star and its luminosity class, but it depends also, 
albeit to a lower degree, on the microturbulent velocity of 
the stellar atmosphere, on metallicity and, within the lumi- 
nosity class of dwarfs, on gravity. A change in Vturb from 
5 to 20 kms ' can shift the spectral type by up to half a 
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Fig. 11. Variation of the nebular fine-structure line strength ra- 
tios [Ariii]/[Arii] 9.0/7.0 and [Nem]/[Neii] 15.5/12.8 yum as 
a consequence of modifying the stellar and nebular metallicity 
from 2 to 1/1OZ0. The adopted SEDs are models #2 (solid tri- 
angles) and #5 (solid diamonds). The dotted line connects the 
models at solar metallicity (star symbols, cf. Fig.|5Jl. A unique 
ionization parameter, U, is adopted for all the models. We stress 
that unless U is constrained, these line ratios cannot be used to 
derive an absolute value for Tgff. 



subclass (to a later type). A change in metallicity from 1/10 
to 2 times the solar metallicity can shift the spectral type up 
to a subclass (to an earlier type). A decrease in logg of 0.3 
dex implies a shift of about one spectral type (to an earlier 
type). 

We confirm the decrease in Tcb for a given spectral type due 
to the inclusion of l ine blanketing. In p articular, compared 
to the calibration of lVacca et al .' ( 1996), which is based on 
plane-parallel, pure hydrogen and helium models, we find a 
decrease which varies between ~ 2500 and ~ 3800 K. The 
decrease is the largest for stars in the spectral range of 05 
to 09. 

The comparison with the stellar features observed in the K- 
band show that variations in the microturbulent velocity do 
not have a significant impact on the equivalent width of the 
Bry, He I and Hen lines. In the case of very hot stars vari- 
ations in the metallicity have an appreciable influence on 
the equivalent width of the He 11 line. In comparing the ob- 
served C rv emission features, important in the spectral type 
calibration using this wavelength regime, the models show 
a discrepant behaviour. C iv is predicted to be in absorbtion 
while it is observed to be in emission. 
The spectral energy distribution between 303 and 450 A is 
best described by the cmfgen SEDs. 

The spectral energy distribution below ~ 450 A is shown to 
be highly dependent on metallicity. This is reflected by the 
behaviour of the nebular fine-structure line ratios such us 
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[Nem]/[Neii] 15.5/12.8 and [Arm]/[Arn] 9.0/7.0 /urn. The 
dependence of these Une ratios on metallicity complicates 
their use as diagnostic tools for the effective temperature 
of the ionizing stars in Hn regions and the age dating of 
starburst regions in galaxies. 
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